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ABSTRACT 


This paper presents a mathematical model of carrier aircraft 
operaticns in anti-air warfare. The mcdel developed herein was a ne- 
cessary outgrowth of a large scale study conducted by Stanford Research 
Institute (SRI) for the Navy. Cne of the problems facing SRI was de- 
termining the number of aircraft stations that a carrier task force 
could fill continuously over some given time pericd. Therefore the 
measure of effectiveness produced by the model is a probability dis- 
tribution over the number of Combat Air Patrol (CAP) or Airborne Early 
Warning (AEW) aircraft able to be maintained on station as a function 
of time in an operating area. 

The medel is developed as a time dependent Markov process with an 
appropriate stochastic matrix of transition probabilities. The main-~ 
tenance process is described mathematically as a pulsed input, multiple 
server queue. The model has been programmed for use on a high speed 
digital computer, the CDC-1604, installed at the U. S. Naval Postgrad- 
uate School, with arbitrarily specified parameter values. The computer 
version of the model with operating instructions is currently available 


at the school. 
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PREFACE 


This paper was started as a summer project by the author at the 
Naval Warfare Research Center, Stanford Research Institute, Menlc Park, 
California in June, 1960. The basic ideas and concepts were formulated 
jointly by Mr. Jack A. Butler, Dr. Donald Guthrie Jr., both of SRI, and 
the author. Further development and translation of the model into a 
computer program was done by the autkor at the U. S. Naval Postgraduate 
School, Monterey, California from July 1960 to March 1961. Although 
the original concept did not include a sensitivity analysis of parameter 
values, the author believes that the development would be incomplete 
without a discussion of input parameters and the effects of varying 
their values. 

A problem similar to the one considered herein has been studied by 
the RAND Corporation, using an Air Force squadron of aircraft with a 
complex maintenance system. The RAND model, reported in RM~2374, is 
much more detailed in maintenance aspects, but it is a Monte Carlo model. 
Consequently a large number of runs with the same set of parameter 
values are required to place confidence in the results. It was our 
belief that the maintenance simulaticn could be Simplified and the en-= 
tire operating cycle be abstracted as a Markov process. In this respect 
we were successful, but there remain questions to be answered. 

Perhaps the most important is "how clcsely does the model approx- 
imate the system it describes?" There is a lack of satisfactory sta- 
tistical data to answer this question, for several reasons: 

1. It is believed that the concept of maintenance "spots" as 

developed in assumption 10, page 6, is a novel but useful one; 
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2. The last full scale operation under combat conditicns 
for which any information is available is the Korean 
War; 

3. The maintenance reporting procedures currently in use by 
the Navy ao not provide the information required for in- 
puts te this study. 

Attempting to answer the above question by collecting and analyzing 

data is considered a (large) separate problem. kHowever, the snswers 
obtained from over 500 runs of the computer model with greatly varied 
inputs always produced answers that were credible. If the collection 

of data from large scale fleet exercises could be undertaken witk this 
model in mind, the results could be compared with operationsel experience 
and an indication of "goodness of fit" would be obtained. Further, it 
might be possible to refine the model to more closely approximate the 
real world system, end provide a dynamic model that could be used by 
operations planners for a wide range of conceivable situations. Even 
in its present state of development, the model seems satisfectory for 
providing quantitative answers to a problem that has previcusly been 
solved only by estimates of the planner. It is believed that the values 
of “upper and lower bounds" on the answers discussed in the sensitivity 
analysis is particularly valuable to the Navy. 

I wish to express my Sincere gratitude to Professor Franklin F. 
Sheehan, who in class introduced me to many of the matrematical con- 
cepts used in this paper, and as faculty advisor provided the en- 
couragement and egsistance necessary to apply these concepts; to Fro= 
fessor Jack A. Borsting for his helpful advice as second reader; to 
Dr. Donald Guthrie Jr. and Mr. Jack A. Butler of SRI, whose names 
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shoulc appear as co-authors of several ideas in this paper; and to my 
wife, who, while enduring many lonely hours of ccmputation and writing, 


provided me with the inspiraticn and clerical assistance necessary to 


do this work. 
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CHAPTER I 


THES FROBLEM 


The broad problem considered in this study is typical cf those fac- 
ing any operations planner. It may be stated as follows: "given a cer- 
tain number and type oi aircraft, end a certain desired coverage of a 
defensive area, what is the optimum stationing policy for these aircraft? 
How many shall be on station at all times, and where should these sta-= 
tions be? If these values sre fixed, what are the chances that the 
coverage will be satisfactory for the period under consideration? What 
are the effects of changing any specific stationing plan? Finally, what 
confidence may be pleced in the answers obtained?" The planner or de- 
cision maker congidergs all information or cate that he considers pertin= 
ent to this problem, and derives answers to the above questions. The 
answers are a functicn cf his experience, estimates, and essumptions. 

If he concludes thet it is unlikely the cesired coverage will be main- 
tained with the given number of aircraft, he may recommend augmenting 
this number to improve the chances of success. 

The model developed in the following pages will provide quantita- 
tive answers to the above questicns. Its development is in the logical 
sequence of: 

1. Consideration of the initial conditicns; 

2. Assumptions; 

3. Inputs derived from the initial conditions end assumptions; 

4. Formulation of the model to translate inputs into outputs; 

5. Sensitivity of the model to changes in the initial conditions 


and essumptions. 
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Although the model was designed as representing aircraft carrier 
operations, it may be readily applied to any type aircraft squadron with 


similar operating procedures. 





CHAPTER II] 


THE MCDEL == INITIAL CCNDITIONS 


The initial conditions presented are those required to charac~ 
terize a typical carrier task force operating at sea either at present 
or some future time. They essentially limit the study to answering 
specific questions as mentioned in the statement of the problem. Con= 
clusions based on the answers obtained are not drawn for any related 
questions, e.g., are more carriers needed in a task force requiring a 
specified defensive coverage?, etc. 

The carrier task force is composed of two CVA-59 class attack car= 
riers and other surface units. A certain percentage of the total deck 
space available on each carrier will be allcted to defensive aircraft. 
Both Airborne Early Warning (AEW) and interceptor type aircraft (CAP) 
will be loaded to conform to this percentage allocation. Standard CVA 
operating procedures will be utilized for flight operations, with ree 
covery shortly following launch of succeeding flights to reduce time 
into the wind as much as is practicable. Flight scheduling will reflect 
the above policy also, consistent with endurance times of specific aire 
craft types. 

A defensive zone will be established for the task force. [In this 
zone, continuous surveillance by AEW aircraft and continuous coverage by 
interceptors is desired. The required aircraft may be provided by both 
carriers jointly, or by alternating "duty carriers". The AEW aircraft 
will be one specific type, as will the interceptors. However, considera= 
tion will be given to using different types of interceptors for the en- 


tire operation, in order tc compare aircraft with different input 
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parameters. 

Task force disposition will not be considered other than fixing 
the general shape of the defensive area in order to determine approxi- 
mate values for the number of aircraft stations to be filled. Similarly, 
surface unit equipment configuration will not affect the stationing of 
aircraft and no attempt will be made to match surface defense capa= 
bility with that of CAP, since this study is cencerned only with the 


air aspects of defense. 
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CHAPTER III 


THE MODEL = ASSUMFTICNS 


In order to assure full CAP and AEW coverage at all times during the 

cycle, aircraft will be relieved on station. 

The percentage of the total CVA=5° deck space allocated to anti-air 

warfare aircraft will be held constant. The quantities of intercep= 

tor and AEW aircraft will be varied within this fixed deck space, 

but the percentage will not vary. The attack aircraft capability 

will therefore be the same regardless of what type interceptors are 

aboard. 

It is assumed that the duration of the CAP and AEW cycles will be 

24—72 continuous hours. The model has been programmed for 72 hours 

to illustrate the outside limit. 

The maintenance crews will have suificient notice pricr to the be- 

ginning of the CAP cycle to complete the 60 and 120 hour checks on 

all interceptor and AEW aircraft. It will not therefore be ne- 

cessary to perform these checks during the CAP cycle. 

It is assumed that no aircraft parts will be deliverec to the car- 

riers during the CAP cycle. This assumpticn is realistic for the 

following reasons: 

ae The task group will be within range of enemy aircraft during the 
CAP cycle which would make delivery very precarious. 

b. Delivery aircraft might assist the enemy in pin-pointing the 
task group's location. 

c. The task group will be far removed from supply points, which 


would make it improbable that parts shortages that develop could 
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be remedied during the 72 hours. 

6. Once an aircraft becomes "Aircraft Out of Commission Parts" (AOCP) 
it will remain in that status for the balance of the CAP cycle 
(see Assumption #5). 

7. Every interceptor on a CAP station will always have its full combat 
fuel reserve so that it can be recruited to intercept raiders if 
necessary. The CAP aircraft will be relieved in sufficient time so 
that they will never need to use the reserve except in combat. 

8. Only minor maintenance will be performed on the flight deck during the 
duration of the CAP cycle. 

9. Except in extreme emergencies, aircraft will be recovered by the 
carrier from which they were launched. Since carriers will be 
loaded with a full complement of aircraft, the recovery of aircraft 
from another carrier would cause severe space problems. 

10.Maintenance for all aircraft of any given type will be conducted as 
follows; 

a. Crews will be formed into a number of groups, and will be as- 
Signed certain areas od the hangar deck in which maintenance may 
be conducted, 

b. Aircraft will be repaired in these maintenance "spots" by these 
crew teams. The teams may be composed of similar or different 
rating groups. 

c. All types of maintenance may be grouped as one class, with no 
distinction among airframes, engines, electronics, or other types. 

ad. Each of the maintenance "spots" will repair aircraft at some mean 
repair rate, say \ , and this mean repair rate is identical for 


each spot. 





e. The maintenance crews will work in these spots for the entire 
duration of the planned 72 hour cycle, probably on a watch basis. 

l1i.The probability that an aircraft will have a failure requiring mainte= 
nance immediately prior to, during, or just after launching (1.e., 
become a "dud" at launch) is identical for similar type aircraft, 
but independent within that type. 

12.The probability that an aircraft will have a failure during flight 
that will require maintenance upon return is also identical for the 
same type and independent within type. Both assumptions 11 and 12 
say that the failure of any specific aircraft will not affect the 
probability that the remainder of its type will fail when launched or 
in flight. 

13.No aircraft will be lost due to enemy action or casualties sboard 
ship. If attrition is expected, the expected number of aircraft 
lost during the operation may be deleted from the initial loading 
and the resulting answers should represent lower bounds on soluticns. 

14.Aircraft becoming AOCP during the operation will be deleted from the 
initial loading before commencement of flight operations, again give 
ing a lower bound. Assumptions 13 and 14 are made to keep the num- 
ber of states under consideration at any given time within some 
reasonable size. If this is not done, the capacity of the computer 
memory for which the mathematical model is programmed will be greatly 
exceeded, and time required for computation would rise by at least 


a factor of four. 





CHAPTER IV 


THE MCDEL == INPUTS (CR PARAMETER VALUES) 


The inputs used in the mcdel were derivec from the initial con- 


ditions and assumptions. For any simulation of real world events, the 


inputs must be realistic and reasonable, or else the development of the 


model is likely to proceed on a weak foundation. Justificaticn is 


supplied where considered necessary in the following discussicn. 


1. Essentially fixed inputs for a given type aircraft: 


Go 


de 


Co 


Time required to climb to altitude and proceed to assigned 
station after launching. 

Time to return from station and be recovered by carrier. 
Endurance time (either on station or total) for a given 
loading. 

Average turn around time to ready for next launching aboard 
carrier. 


Percentage of deck space occupied by each aircraft. 


2. Inputs that must be estimated: 


Ge 


b. 


The probability that any given aircraft will incur equip= 
ment failure while airborne requiring non-routine maintenance 
upon return to carrier. 

The probability that any given aircraft will experience 
equipment failure just before, during, or after launching 
that will cause the flight to abort, and the aircraft to 
undergo non-routine maintenance. (The probability of a 

"deck dud"), 


The number of maintenance "spots" available for repairing 


8 





aircraft on the hangar deck. 

d. The average time to repair failed aircraft and return them to 
a fully ready status. This may be given as a mean repair rate. 

e. Number of aircraft becoming ACCP during the operaticn. 

f. Percentage of aircraft complement available at beginning of 
Cperation in a fully ready status. 

3. Arbitrary inputs fcr a given type aircrafts: 

a. The number of stations desired filled continuously during the 
operation, 

b. Intra-carrier operating doctrine, i.e., alternating launches 
by each CVA or simultaneous launches of approximately half the 
required number. 

The fixed inputs require only minor discussion, since they may be 
obtained from aircraft performance and characteristics charts, and from 
deck space allocations mentioned in the assumptions. If the aircraft 
is still in the design or development stage, planned values may be used 
initially and revised when actual performance data are available. 

The estimated inputs must be fixed or some realistic range of values 
given. For operational aircraft, statistical information might be used 
in eStimating the inputs, but for future types: or models, cnly estimates 
and drawing beard specifications are available. Stanford Research In- 
stitute has done considerable work in attempting to determine a range of 
values for these inputs for specific types of aircraft, using manufactur= 
er's specifications, Bureau of Weapons estimates, and aircraft squadron 
operational experience. The sensitivity analysis discussed later in this 
paper will cover in more detail the effects of varying these inputs. 


It is the author's belief that a series of computer runs must be made 
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for each set of fixed and arbitrary inputs, covering a range of values 


for the estimated inputs. 
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CHAPTER V 


THE MCDEL -- DESCRIPTICN 


The Air Cperations model described herein takes certain predeter- 
mined inputs and operates on them to produce one output which is a 
direct measure of effectiveness of the air operations system for one 
particular type of aircraft. The output is presented as a probability 
distribution over the number of aircraft that can be maintained on 
station for a given length of time, given these inputs. Although the 
model is probabilistic in nature, it is in no way random or "Monte 
Carlo": it produces an analytical solution for the probability distri-~ 
bution. Therefore, if reasonable confidence can be placed in the in- 
puts, the results obtained will reflect the accuracy of the model in 
simulation of the operations system in question. 

Basically, the mathematical model is formulated as a specific type 
of stochastic process called a Markov chain (2 |[3|[4 J e suppose we 
examine at specific points in time, t, a process that proceeds contin-~ 
uously in time. Let the outcome (state) at each observation time be one 
of a given fixed set of outcomes B) (k= 1,2, ...n). Let us assume 
that when we started, there was a probability distribution over the 
possible states, called alto? S [a §*?, a, te), a gto) |, SO 
that at the start the probability of being in state KE, will be q,, %0) ; 

Let us define a conditional probability Pi3 as the probability 
that at observation time t+1 outcome é, was observed, given that at 
observation time t, outcome E; was observed. If we arrange these con- 
ditional probabilities in a square matrix P called a matrix of trans- 
ition probabilities, it will have the following properties: 
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1. The elements are non-negative (p3 420) for all i and j. 


T\ 

2. The row sums will be unity ( jp 
jel 

A matrix satisfying the above conditions is called a stochastic or 


», poe )iforvalind. 
iJ 


Markov matrix. A Markov chain or process is completely defined by a 
Markov matrix of transition probabilities P and a column vector gto) = 
[a, to) qo '*o), ae an '*? | which gives the probability distribution 
over the states at the start. Note that the dependence of the Markov 
Chain extends over one observation period only. 

To establish the states B, for the model, let us consider two 
random variables defined as follows: 

a. X,(t) = the number of a/c flying at time t, not having flown 

in - previous launch-to-launch interval. 

be X(t) = the number of a/c in or awaiting maintenance at time t. 
Then if we consider the vector EXGr Xp(t) | as a pair of random 
variables, we will have a bivariate stochastic process with the possible 
states E. characterized by pairs of values [2]. If we have N total 
aircraft of a given type aboard the carrier, and we wish to maintain 
A aircraft on station, the range of values of BE, will be from (0,0) 
to (A,N) or Be (xsy)s O€xSA, OF yEN. 


The initial probability distribution alto) over the states E, is 
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simple to describe, but determination of the specific distribution may 
not be easy. Since it is assumed that a task force will anticipate im- 
pending operations and enter an objective area with as few aircraft in 
maintenance as possible, the most likely value for X(t) isezeres Te is 
also assumed that aircraft operations will not commence until entry into 
the objective area, at which time full cycling will start. Therefore the 
initial distribution g{to) is assumed to be q, to) 2 one for & = (0,0) 
and q, ‘to = zero otherwise, but the model has provisions for considering 
any other possible initial conditions at the discretion of the operations 
planner. 

Let us define a typical a/c operating cycle as a unit interval of 
time. Cbservations of the process sare made at the initial launch and at 
successive unit intervals of time (launch times). To develop the con- 


ditional probability, » Of the system going from state E; to state Bs 


a i 
in a unit interval of time, let t be a typical launch time. Launching 
Operations commence at time t and continue to time t + e. At time t+ T, 
aircraft from the previous launch are recovered. At time t + 1, the 
succeeding group of aircraft is launched, and the cycle repeats. ~ 
Therefore, the unit of time is the launch-to=-launch interval, t to t+ 1. 

In the interval t to t + e, aircraft are launched with some 
probability of success. Those that are not successfully launched enter 
maintenance, In the model, this launching process continues until the 
desired number of aircraft (A) are successfully launched, or until there 
are no more available aircraft on deck. Although this time interval is 
of finite length it is small compared to the unit of time and it is 
treated as a point in time by letting e become zero. 


Commencing at time t + T, the aircraft from the previous flight are 
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recovered. Again this interval of recovery fromt-+Ttot-+-T+e is 
idealized as a point in time. The consequences of assuming that the 
launch and recovery intervals are of zero time duration will be discussed 
in more detail in the sensitivity analysis. When the aircraft are re- 
covered, some may have experienced equipment failure in flight charac= 
terized by a failure probability described in assumption 12, and will 
also enter maintenance. Those that do not have failures will be rearmed, 
refueled, and be made ready for succeeding launches. 

During the unit time interval, maintenance will be performed on 
those aircraft in a non-eready status, and a certain number of aircraft 
will be repaired as noted in assumption 10, There is a chance that the 
entire unit time interval will not be available for maintenance, This 
possibility is discussed in the sensitivity analysis. 


a(to) 


To summarize, we start with some initial distribution, 5» over 
the possible states B, at time t,; launch, recover, and repair aircraft 
in a unit time interval; then repeat the process for the succeeding in- 
tervals until the end of the operating psriod (see figure 1, pg. 16). 

if we know the transition probabilities for all states of the system, we 
may form the transition matrix P, the elements P(«,i), (B55) being the 
probabilities of going from state («,i) to state (6,3) in unit time for 


all possible states, and thus determine the probability distribution over 


the states (853) 2 at the start of any succeeding cycle as follows: 


t+1 a + 
a! = Pr 2 (t +1) Ss (@53) |= a ) ° Fp ° 
Finally, the desired measure of effectiveness may be cbtained at any 
unit time t (i.e., the start of any cycle) by summing appropriate main- 


tenance state probabilities, i.e., the probability that « a/c are in the 


14 





7 


N 
air: Pr[x,(t)<2«]e 2, 
=O 


Pr [X(t) = («,4) |. 


td 


This is the probability that « aircraft are launched regardless of the 


number in maintenance. 
The transition matrix P 1s complex in nature and the development 


of its elements is described in appendix A. 
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Basic Diagram of the Cperating Cycle 
Figure 1 
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CHAPTER VI 


ij 
RESUME OF COMPUTER PROGRAM 


The computer program of this mathematical model was written for the 
CDC = 1604 digital computer. The following items are available at the 
U.S. Naval Postgraduate School, Monterey, California: 

a) IBM Cards of the program in CDC CO = OP Assembly Routine (AR) 

language, 

b) A complete line printer printout of the assembled program. 

c) A punched paper tape of the program in "AR" format, ready for 

assembly. 

ad) Description of the program and Operating Instructions. 

e) Flow charts of the program. 

f) Sample copy of output for a specific set of input parameters, 

g) The complete program in Binary Coded Decimal format on magnetic 

tape. | 
Readers interested in obtaining any of the above items or information 
concerning them should address inquiries to Secretary, Department of 
Mathematics and Mechanics, U.S. Naval Postgraduate School, Monterey, 
California. 

The program was written in several parts to facilitate checking out 
and "debugging". Initially, the input parameters for a specific run or 
series of runs are placed in a table in memory. The required values of 
Ps; (T) are then computed in two separate subroutines, one for direct 
evaluation of equations 3 and 4, pg. 29, and the other for solution of 
the system of differential equations (equation 1, pg.32) using numer-= 


ical integration. The values obtained are stored in a table in com= 
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puter memory, using a perallel table storage of the values (i,j) to 
facilitate retrieval and retain the accuracy of the computed values. 
Next, the values of Maan) and Yepn are computed, using subroutines 
included in the program, and stored in memory using the same procedure 
described above for Pi (1). The program then generates the transpose 
of the transition matrix, P'y and stores the elements in memory in cone 
Secutive order, The matrix is generated in transposed form to facili- 
tate indexing and future matrix multiplication. A jump is then made 
to a routine that sums each row of the matrix and transfers these sums 
in floating point decimal format to magnetic tape for printout to 
check the accuracy of computations. 


(t,) 


The initial probability distribution vector Q is then formed, 
the matrix product gto} o PF x gto +1) is generated, and the measure 
of effectiveness is obtained from this new vector (by summing over all 
values of i for fixed oe ) and transferred to magnetic tape in floating 
point decimal format for subsequent printout. This new vector alto 1) 
then replaces the original vector gfto) and the next matrix product is 
formed. This process continues wuntii the end of the operating period 
is reached; then if investigation of other initial distributions is 
desired, the above process is repeated. If not, computation of the 
measure of effectiveness is completed, and (at the discretion of the 
operator) either the program stops, or the last value of the measure of 
effectiveness Pr x, (t) so al is checked to see whether it is less 
than some preset value (presently .9). If not, A is increased by one, 
and computation is resumed, omitting certain parts of the program that 
do not require re=running. Finally, when the above condition is satis- 


fied, computation stops. 
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There is a separate section of the program that may be used if it 
is desired to conduct a sensitivity analysis similar to that described in 
this thesis. Initial values are inserted and computation is begun in 
a different location. Parameters are varied within the same range of 
values shown in Appendix C until the series of runs is completed. Ree 
sults are printed in consecutive order on magnetic tape for printout 
on a line printer. 

The computer program with all subroutines and storage cells except 
those for decimal conversion and output is approximately 3000 two in= 
struction words in length. Essentially all computations sre done using 
floating point arithmetic. There are 28,160 cells (of a total of 32,768) 
reserved for matrix elements, restricting the values of A and N in that 
(A £ 1)(N 4 1) £ (28,160)® xe 167. If the values of A and N should not 
Satisfy this restriction, it would be necessary to re-program the model, 
perhaps storing portions of ths matrix on magnetic tape while generat-= 
ing the remainder. 

The size of the matrix is the major factor determining running time 
for the program. For an A of two and an N of five, running time is 
approximately two minutes, while for an A of six and an N of 20, this 
time increases to approximately 35 minutes. Minor attempts to rewrite 
portions of the program and reduce the length of time required for com= 
putation have not significantly lowered these values. Although the 
author is not an experienced programmer, it appears that the number of 
arithmetic operations required for generating each element of the matrix 
cannot be reduced, and these approximate times must be accepted. Howe 
ever, it is emphasized that the cost of 35 minutes of computer time 


compares favorably with the cost of 72 hours of task force time. 
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CEAPTER VII 


SENSITIVITY ANALYSIS 


The model that has been created to transform inputs into outputs 
does not appear susceptible to a mathematical sensitivity analysis due 
to its complexity. However, a lengthy program of variation of input 
parameters was conducted to compare numerically the resulting variations 
in output. To limit the analysis to inputs which were essentially es= 
timates, the following values were fixed for a simmlated type aircraft: 

As Four a/c desired on station a/c = aircraft 

Ns 13 a/c available aboard CVA 

Dz: Six maintenance spots available 

Cycle length: Frur hours 

7: 02 of cycle length 

Duration of operations 72 hours 

The following inputs were varied within the ranges showns 

The probability of equipment failure in flight (p): 0(.05).5 or 

O(.1).5~ 

The probability of a deck dud (P,): 0(.03).21 or .0(06).18. 

The mean repair rate (4): 4(4)20. 

The fraction of cycle time available for maintenance 1.0(.2).8. 

Tabulated and graphic results of these runs are shown in Appendix C 
and D and are to a large extent self explanatory. They illustrate how 
variation of each parameter affects the measure of effectiveness. 

Although inputs may be varied at the discretion of the user, it is 
envisioned that basically two types of group runs of the computer model 


would be made. They are: 





1. Predetermined fixing of the probabilities and maintenance 

values and variations of A, N, and Cycle Length (and with it 
TY), to determine the optimum loading, stationing, and cy= 

Cling of a particular type of aircraft. 

20 Fixing the number of aircraft available, their stationing 
and cycling, and varying the maintenance values and the prob= 
abilities of failure, as has been done for the results pre- 
sented herein. 

The first type of group runs was conducted for SRI using current 
and future fleet aircraft parameters in order to determine optimum 
loading stationing and cycling policies for these types of aircraft. 
The results cannot be included due to their classification, but real- 
istic output values were obtained. 

For any mathematical abstraction of a real world process, the 
question of goodness of fit of the model to the actual situation is the 
basic test of the value of the model. As mentioned in the introduction, 
this question has not been answered directly; however, the answers ob= 
tained for any given situation "seemed reasonable", One of the funce= 
tions of an operations analyst is to recommend or arrange for suffi- 
cient collection of data, usuaily to provide quantitative values for a 
statistical analysis of a problem. Therefore, it is the author's re= 
commendation that if this model is used by the Navy (and he hopes that 
it will be), a program for collection of required data be established 
for future large-scale carrier operations. Perhaps some of the required 
inputs are available from already existent reports, and the remainder 
might be obtained by modifying their format. A proposed data collection 
form for inputs required by this model is given in figure 3. 
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MCNTH _ YEAR _ SQUADRON __sss—té‘“éT'YPES:«OAA/C 


Day 





1. # a/e aboard 











2. # a/c in ready status 








Sau a/c ACCF 


os Soe mae coc 





4. # sorties 








5. Average duration of 
each flight 





6. # of aborted flights, 
or a/c downed just 
prior to or during 
launch 





as = oer a a 








7. # a/c down for 
routine maintenance 
or modifications 





8. # a/e entering non- 
routine maintenance 
this day 

a. upon return 
from flight _ 

bo from a ready 
status 





CR COO ee 

9. # a/c leaving non= 
routine maintenance 
this day 





10. Time required to 
repair each of these 
a/e (list time in hrs. 
for each a/c) 






ee ne <n S pe Se 


11. # of man hours required | 
for repair (list for 
__each a/c) 


Ee ee ee ee 


Sample Maintenance Report Form 
Figure 3 
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Variation of the length of time available in any one cycle for 
maintenance was considersd for two reasons. First, it appeared likely, 
although not mandatory, that maintenance work might stop during launche 
ing and recovery operations. Secondly, the assumption that launch and 
recovery intervals could be idealized as pcints in time would have to be 
modified if maintenance work was halted in these intervals. For the 
four hour cycie length chosen for the sensitivity runs, a significant 
degradation of maintenance time of 48 mimates was chosen to show the 
effects of continuous as opposed to interrupted work. If the carrier 
were operating only this or a similar type of a/c and maintained four 
hour launch intervals, 48 minutes appears excessive. However, if other 
a/c types of limited endurance were being launched and recovered, this 
value is esaneea to about 15 minutes per launch-recovery operation. 

Since maintenance work might also be interrupted for relocation of a/c 

on the hangar deck, it seemed desirable to consider this longer time. 

On the other hand, for a cycle length of 90 minutes, this degradation ree 
sults in a loss of 18 minutes of time for maintenance work, again a 
likely value. 

The effects of this loss in repair time may be noted in figures 17 
and 18, Appendix D. As shown, the effect of not having the entire in- 
terval available is essentially negligible for a fairly short mean repair 
rate and Pelisr failure probabilities. As the average length of time 
required to repair each a/c increases, the fraction of time the maintez- 
ance system is not working becomes more significant, certainly an ex= 
pected result. The same effect is evident when the failure probabilities 
are increased. 


Fortunately for a mathematical model of this size and complexity, 
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there are certain built-in checks of the correctness of the computed 
answers. Although computational subroutines were checked before wse to 
ensure their exactness, the row sums of the transition matrix P give an 
adequate check of the correctness of a major portion of the program. At 
no tims during the computation did these row sums differ from one by 


Ree and in most cases this difference was less than 107°, 


more than 10 
The small errors, when appearing, were traced to the numerical integra— 
tion subroutine for Pi3h))s and appeared to be the result of inherent 
roundoff errors in the computer, when converting fractions to fixed 
length octal numbers. Considering the number of arithmetic operations 
performed in the program, the exactness is certainly satisfactory. A 
second check is obtainable from the measure of effectiveness, in that the 
sum of the probabilities of launching zero up to A a/c should be one. 
Again, the same exactness indicated above was obtained. Note that the 
computer program produces "exact" answers that are arithmetically cor- 
rect, as opposed to "“accurate™ answers, which to the author imply 


correctness in the sense of duplication of the real world system. This 


accuracy has been discussed elsewhere in this paper. 
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APPENDIX A 


DEVELOPMENT OF THE TRANSITION MATRIX P 


The elements of the transition matrix, Pyyo are the probabilities 
of transition from state Be to state BE; in the unit time interval. Gon= 
Sidering a typical unit interval, at the beginning of the interval there 
is a probability distribution over all the possible states. Let us deo 
fine Yoh as the probability that of f aircraft on a carrier deck ready 
to fly, g will be launched and h will be duds requiring repair, each a/c 
having the same chance of successful launch characterized by a success 
probability s, which is one minus the probability of a "deck dud", Py ° 
Since each a/c launch is in itself an independent event, and we desire 


to obtain A successful launches, ¥¢o4 will have the following valuess 


0 if g>A 

2. 0 ae | 1 BE: 

3. 0 , ifgéaéA, gr+hef 
i 4 = 

Lo (7 3 @! _ s)t 8 if g ZA, g + h = f 
ree, bh ¥ 7 

56 Oe } 3& (1 = g) if geq=A. | | 


The reasoning for these values is as follows: 
1. It is desired to launch no more than A aircraft. 
2. It is impossible to launch and send into maintenance more a/< 
than the number available. 
3. If the required number A are not obtained and there are still 
ready a/c left, launching will continue. 
4. The standard binomial "success-faiiure" expression for indepen— 


dent trials, where all ready a/c are used, but the desired 


at 


number A are not launched. 

5. The standard negative binomial expression for obtaining eg 
successes in g+h- 1 trials, i.e€., the desired number are 
iaunched without utilizing all available ready aircraft. 

Considering next the time t + T when say « a/c return and are reo 

covered, we must determine the prokability that say, m of these a/e 
require repair upon landing and enter maintenance, called We (m). Since 
@ach flight is an indevendent event as far as failures are concerned, 


VWy.(m) will be a standard binomial distribution, i.€., 


VW. (m) = (=) ei = pies m > where p is the probability of equip= 
ment failure in flight. Those a/c that do not fail are refueled, rearmed, 
and prepared for succeeding flights. 

Turning to maintenance procedures, let us define P, (7) as the 
probability of repairing (i = j) a/c in time , given i a/c in or 
awaiting maintenance at time t. An alternate definition would be the 
probability of having j a/c in or awaiting maintenance at time t+, 
given i at time t. Under the assumptions given in part 10, a pulsed 
input, exponential holding time, multiple server queue [5 | may be devel= 
oped as followss let D be the number of maintenance "spots" available, 
These are service channels in queue terminology. Each of these service 
Channels will have an identical mean service rate, ), and an exponent 
ial service or holding time. OEG Report 585 (Confidential) verifies the 
exponential service time for maintenance of aircraft aboard carriers 
during the Korean War, the latest information available for large scale 
Operations. 


Then for each channel occupied, the probability that it remains 


*y 
< 





cS iS 
occupied throughout a length of time IY is e >» and the probability 
that it becomes free (repairs the a/c) in this time is 1 - en AT 
The queue equations which were arrived at independently are 


identical to those given in 6 | o They ares 


dP 
i fon +. oAP, p(t) ~ (n+ LAP, na pit) for O4neD, 
| AP | t ss 
ee Baga !t) “ DAP, ,(t) = D dP, a» yt) for neD, 
dt ‘ 
with P, 460) =. one if i= j, zero otherwise. If there are fewer air 


craft awaiting or in maintenance than there are "spots", there is no 
queues, and the resulting distribution of Py? is binomial, since each 


Spot is considered to work independently of the others. In this case, 


. HAT tes 7 ped ' 
36 P(e) 2 fl) (laoae AT): a ee 4 e 
J J 


Naturally, if i<j, Pk?) is zero since thers is no way for the pulsed 
queue to increase with no input in the interval of repair considered. 


Equation 2 above may be solved in clossd form to give: 
) 


leon a 
fs re (t) = (DAt)? = | e DAt for! i2 n 2D. 


25D (i = n)! 

The solution of the system of equations given in 1 above is dis= 
cussed in Appendix B. Since the work involved in producing and eval- 
uating a closed form solution is at least comparable in time to that 
required for solution by numerical integration, a program for solution 
using numerical integration by the CDC=1604 computer was written and 
used, and provided completely satisfactory rapid answers. 


To summarize the values of Pi 3(™) then, : 


oa 





0 i oo a 
(=) ae er AMI =» j sia 246 1a 


P'* pated par . > of 
pe oa if "iy ope Vee 
Numerical solution of equation 1 ioe Jie sd, a 
To generate the transition matrix P, we must consider all the 
possible events taking place in the unit time interval as follows: 
At time t, we start with « a/c in the air and i a/e in maintenance; a/c 
are launched and enter maintenance at time t + e3 a/c are repaired from 
time t + e until time t + T, when a/c are recovered and some may enter 
maintenance; a/c are repaired in maintenance until the end of the unit 
cycle, t+ 1. Then in words, the probability that there are g a/e fly- 
ing and j a/c in maintenance at t+ 1, given «a/c flying and i in mainte 
enance at t is the probability that the system changes from state («,i) 
to state(@,j) in unit time. Each Pay is composed of the triple sum 


over all possible values of the product of four event probabilities, i.e. 


P(,4), (855) se Vie (#4) 651 ; Pi41,\(T) Tem) : Pom, 32 = T). 
. )™ 
Explanation of the indices is as follows: 


ee Vive (ocei) 581 
which 6 are launched, and any number, 1, from 0 to N=(K+ i) “G o@ 


» There are Ne(«#+i) a/c available at time t, of 


enter maintenance. 
b. Pied yi T)> Initially, there were i a/c in maintenance, to which 
d 
1 a/c are added when launching occurs. In time T, the maintenance 


queue diminishes in length to k as a/c are repaired. 


Ce We (m). At time T, m of the « a/e flying enter maintenance, 


30 





where m may range from 0 to &. 
P (i =T). At time T, the maintenance queue recsives m 
ke mj 
additional a/s, and in the remainder of the interval, (1 ~ T), 
the queue diminishes in length to j s/c. The limits on k are 
from zero to the total number of a/c that could enter mainten- 
a 
ance in the intsarval considered, i + l. 

The transition matrix is of order (A+1)(N*1) x (A+1)(N+1) 
since all possible combinations of («,i) are allowed. However, ml 
states occur in the following cases: 

1, H+ i PN 3 

ame + | pall a 

Zo K+14@B > No 


In each of these cases, it is impossible to use mors aircraft than there 


are available, and consequently all these probabilities are zero. 


ek 





APPENDIX B 


SCLUTICN CF TEE QUEUEING EQUATICNS 


The equations of balance for the pulsed input queue (equations 1 


and 2, Appendix A) are 


UG 


f 1. oF sa) o nAP. a) (n+ LA, « Wee) forsOmm< D, 


dt +9 
2, —ignt®) r A | 
3 + -- py, Ps nt) = D Pi one tt) for n2D, 


with Br uy s= one if 1- j, zero otherwise. As noted in Appendix A, 
equation 2 is solwable in closed form and this solution is given there. 


When n= D-= 1, equation 1 above becomes 


dhe Bt (4) > 
° _ i gDe— ae => P D> ec = Pp. a © 
at pa i,por®) = °4 i,pft) 
, i=D .-DAt 
From Appendix A, P, .(t) = (DAt)? e , 
pp 9 i,ps Js (pat dn. 
Denoting P. (t) as x, equation 3 is of the form s% + Qx «= Ry, 
de 


where Q and R are functions of t only or constants: This is a linear 


differential equation and its solution for x is 


xe R 6 dt +6, © a constant of integration +t) 


(a at ( fQ dt 
| 


Satisfy the initial or boindary conditions. Substituting for Q and R 


and simplifying, this becomes 


fe ot Oand) +5 P pag) cet Di) tere 
— a | 


which eventually reduces to 





2 " 7 =. [" A a 
er DAt pe D+ 1 | ° Se ke y ( t)J | as a solution 
me! 





for P, D z(t). Note that the original limita of i and j for Pi yb? 
9 7 F 
were j<D4i for this system, so that the first value of i, D+rl, gives 


4a solution of 


SpAt” eaeP Re 
D re »~ l= At dl To give a numerical example, 


Xe e 
suppose D » four maintenance spots, A is one half the cycle length in 
time, and t is one fourth the cycle length in times then the probabil= 
ity that the maintenance queue would change from five (D+1) to three 


(D-1) a/c awaiting or in maintenance in one fourth the cycle is 


a 
P, (4) = eo uy? LeV®.1-1/8 | = .079 to slide rule 
9 


accuracy. 


For the next smaller value, n= D= 2, equation 1 becomes 


aP. (t) 
~AsD -2 ¢ (D= 2)AP, p(t) = (D=I)APy p _ y(t). 


Using the expression for P, D 1 {t) given above, the solution for 


Pip _ git) is 


| v-D , 
Pip _ 3(t) x en (D =2)At A(D=1) pa fe . At -F ae) | 


mJ 
dt + C6. 

This integration has been carried out and an expression for Pip - g(t) 
obtained in explicit form, as has been done for Psp A 3 ft) Both of 
these expressions are lengthy and complicated, and succeeding expressions 
for decreasing n will be even more complex. Since the expressions must 
be coded for computer solution for utiligation with the remainder of the 
model, it appeared that solution of the system of differential equations 


could best be handled by direct numerical integration using a subroutine 


of the Runge-Kutta-Gill type available for the CDC-1604 computer. The 


oo 





resulting program provided rapidly obtained accurate values for Py y(t) 
for all requisite values of i, n, and t, using the initial values of P 
and P' at time zero and a delta t of .001 or .005. In addition, it is 
interesting to note that corresponding values tabled in K: | were 


verified. Accuracy of the solutions obtained is discussed briefly else 


Where in this paper. 
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TABLES OF NUMERICAL RESULTS OF CCMPUTER RUNS FOR VARICUS INPUT PARAMETERS 


APPENDIX CG 


TABLE 1 


1. Mean Repair Rate A equals four hours. 


Ao 


PY 


20 
065 
o7 
015 
of 
085 
oJ 
o95 
1.0 


NOTE: 





be 
Pr 

200 

ol 

02 

o3 

oh 


Entire interval available for maintenance. 


-00 
©9953 
9927 
9893 
09847 
9788 
9714 
» 9622 
9511 
agate 


20 
09902 
09862 
o9810 
9745 
9665 
569 
9454 
0 9320 
o9166 


206 
09822 
9764, 
9693 
»9606 
09504 
09385 
09248 
09092 
» 8916 


oO 
09705 
09626 
o9oee 
09425 
09301 
09160 
»9003 
o OG 
S050 


olla 

09545 
oSAAL 
9330 
09200 
09056 
» 8896 
8719 
08531 


»8320 


Above runs were not conducted for p 


od 
9338 
09216 
09080 
8932 
0 8769 
0 8592 
» 8402 
o 8195 
0 IYTL 


oL8 
o90E2 
» 8939 
o8786 
Geel 
S443 
08254 
» 8050 
o FOpe 
° 7605 


o55 since there is 


Significant change in the measure of effectiveness. 


06 of interval available for maintenance. 


1.0 
1,0 
LcOmm 
s7792 
0 9980 





1.0 

so 
09986 
09959 
9891 


Or 





m2 


29988 
9961 
09906 
29807 
09649 


9901 
o 9805 


09659 
09459 
09184 


gee 

8774 
» 8616 
S448 
o 8270 
» 8082 
o 7881 
» 7670 
0 14,50 
o 7216 


not a 


— 








TABLE 2 


Mean Repair Rate \ equals eight hours. 
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TABLE 3 


3. Mean Repair Rate Aequals 12 hours. 


a, Entire interval available for maintenance. 
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4... Mean Repair Rate A equals 16 hours. 
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TABLE 4 


Entire interval available for maintenance. 
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TABLE 5 


5. Mean Repair Rate A equals 20 hours. 


a. Entire interval available for maintenance. 
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Be o8 of interval available for maintenance. 
These runs were not conducted since their inclusion would not add 


Significantly to the results presented herein. 
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